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Abstract 

The metric projection onto an order nonnegative cone from the metric projection 
onto the corresponding order cone is derived. Particularly, we can use Pool Adja- 
cent Violators-type algorithms developed for projecting onto the monotone cone for 
projecting onto the monotone nonnegative cone too. 

1 Introduction 

The metric projection onto convex cones is an important tool in solving problems in 
metric geometry, statistics, image reconstruction etc. In almost all applications the pro- 
jection onto a convex cone is part of an iterative process, hence its efficiency is of crucial 
importance. 

Projecting onto an order cone jl] is a fundamental tool for solving isotonic regression 
problems (see [7]). We will call an order nonnegative cone the intersection of an order 
cone with the nonnegative orthant. An order nonnegative cone correspond to an isotonic 
nonnegative regression problem which we will define from an isotonic regression problem 
by superimposing the nonnegativity of variables. 
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An order nonnegative cone is a pointed cone which is the subcone of the corresponding 
order cone which is not pointed. 

A special case of the isotone regression problem is the case of the regression with 
respect to a complete order ( [7], PQ). The corresponding cone is called monotone cone. It 
turns out that the corresponding monotone nonnegative cone is important in the metric 
geometry and the image reconstruction. Whereby the importance of projecting onto this 
cone. 

A simple finite method of projection onto the so called isotone projection cones (cones 
having the property that the metric projection onto them is isotone with respect to the 
order relation defined by these cones) proposed by us (see [6]) has become important 
in the effective handling of all the problems involving projection onto these cones. The 
monotone nonnegative convex cone used in the Euclidean distance geometry (see [2]) is 
an isotone projection one. Our method has become important in the effective handling 
of the problem of map- making from relative distance information e.g., stellar cartography 
(see 

www . convexoptimization . com/wikimization/index . php/Pro j ection_on_Polyhedral_Cone 

and Section 5.13.2.4 in 0). 

Due to the importance for the regression theory of the projection onto the monotone 
cone, there are various efficient methods of projection onto it. These methods emerge 
from the so called Pool Adjacent Violators (PAV) algorithm, which has nothing to do 
with the geometric approach of [6j. The PAV algorithm exploits the specific feature of 
the monotone cone, and due to its efficiency it would be desirable its adaptation for 
the monotone nonnegative cone too. Our note aims to do this by joining PAV and the 
geometric approach specific for [6]. 

In this note we show that the projection of a point onto an order nonnegative cone is 
the positive part (with respect to the lattice structure of the nonnegative orthant in M. m ) 
of the projection of the point onto the corresponding order cone. 

2 Preliminaries 

Let W be a convex cone in M m , i. e., a nonempty set with (i) W + W C W and (ii) 

tw c W, V* g JR+ = [o,+oo). 

The cone K C R m is said pointed, if K n (-K) = {0}. 
The polar of the convex cone W is the set 

W x := {y G R m : (x, y) < 0, V x G K}, 

where (•, •) is a scalar product in IR m . 

If W is a closed convex cone, then from the extended Farkas lemma (or bipolar theo- 
rem, see e.g. Theorem 14.1 in [5] p. 121) is 

W ±L = (W L ) L = w. 
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If Z is another closed convex cone, then W and Z are called mutually polar if Z = W 
(and hence W = Z L by the lemma of Farkas). 

The scalar product (•, •) defines a metric d on W 71 by setting d(x, y) = (x — y,x — y} 1 ^ 2 
for any x, y G M m . Denote by Pw '■ R m — >■ the projection onto the closed convex cone 
W (or the nearest point mapping), which associates to x G M. m its (unique with respect 
to the metric defined by the scalar product (•, •) [9]) nearest point Py/X G W in W . 

The projection mapping P^y onto W is characterized by the following theorem of 
Moreau [5]. 

Theorem [Moreau] Let W, Z C M m 6e two mutually polar convex cones in H. Then, 
the following statements are equivalent: 

(i) z = x + y, xGW 7 , y G Z and (x, y) = 0, 

(ii) x = Py/z and y = Pzz. 



3 Order cones and order nonnegative cones 

Suppose that R m is endowed with a Cartesian coordinate system, and x G M m , x = 
), where x % are the coordinates of x with respect to this reference system. 
Endow the index set {1, . . . , m} with a partial order ^. The order cone [1] with respect 
to the partial order ^ is defined by 

={iG M. m : x l < x j whenever i ^ j}. 
The order nonnegative cone corresponding to the order cone is defined by 

= n R™. 

We obviously have 

if^C% (1) 

4 Projection onto an order nonnegative cone via the 
projection onto the corresponding order cone 

Let b = (b l , . . . , b m ) G M m and w = (u>i, . . . , w m ) G M+ + a vector of positive weights. 
The partial order ^ defines a directed acyclic graph over the nodes {1, . . . , m} such that 
G whenever i ^ j, where E^ is the set of edges of the graph. Recall that the 
isotonic regression problem is the following minimization problem: 

Minimize Ya=i w i( x i - h) 2 
subject to Xi < xj, G E^. 
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We define the isotonic nonnegative regression problem as the following minimization prob- 
lem: 

Minimize Yn=i w i( x i ~ h) 2 
subject to Xi < Xj, V(z,j) G E-<, 

Xk > for k = 1, . . . , m. 

If we define the scalar product (■, ■) by (x, y) = wiXiyi + ■ • • + w m x m y m for any x, y G R m , 
then the isotonic regression problem is equivalent to projecting b onto W-<. Similarly, the 
isotonic nonnegative regression problem is equivalent to projecting b onto K^. However, 
the following results hold for any scalar product (-, •) in R m . 

For any vector a = (a 1 , . . . , a m ) G R m denote by a + = sup{a, 0} its positive part with 
respect to the lattice structure of R™. Thus, if m = 1, then a + = max{a, 0} and if m > 1, 
then a + = (a 1+ , . . . , a m+ ), where a l+ = {a l ) + for any i G {1, . . . , m}. Similarly denote by 
a~ = sup{— a, 0} its negative part with respect to the lattice structure of R™. Thus, if 
m=l, then a~ = max{— a, 0} and if m > 1, then a~ = (a 1- , . . . , a m ~), where a l ~ = (a l )~ 
for any i G {1, . . . , m}. It is easy to see that a = a + — a~ . 

Using the notation in the preceding section we have the following result: 

Theorem 1 For any y G R m we have 

P K ^y = (Pw^y) + - 

Proof. We have from Moreau's theorem that 

y = Pw^y + P w ^y^ ( p w^y, p w ^v) = °- ( 2 ) 

Denote v = (Pw^y) + ■ Using the notations in the theorem, let u = Pw^y — v, that is, 

u = -(P w ^y)~. 
Denote z = Pyy^V- Then, fl2]) becomes 

y = u + v + z with (u + v, z) = 0. (3) 

First we show that from the special form of W^, we have it, v G W^. Indeed, if we 
denote x = Pw^y, then 

Uj — X j ... j X ^ 

and 

v = (x l+ ,...,x m+ ). 

It is easy to see that the functions R 3 t i— > —t~ and R 3 t i— > t + are monotone increasing. 
Since x G W-t we have x l < x^ whenever i ^ j. Hence, by using the monotonicity of 
R 3 t i— > —t~, we get u l < v? whenever i -< j. Similarly, by using the monotonicity of 
R 3 t i— > t + , we get v l < u 3 ' whenever i ^ j. Thus, u, v G which together with ([3]) and 
z G yield 

( u ,z) = (v,z) = 0. (4) 
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From C (see (JTJ)), it follows that 

W$ C K^. (5) 

From the fact that all the coordinates of u are nonpositive and the elements in 
have nonnegative coordinates, it follows that 

u G K$. (6) 

Let us write now 

y = v + (u + z) 

and observe that v G K-<. Further u + z G K 1 - from ([5]) and (JBJ). We also have 

(t>, it + = (f , u) + (v, 2) = 0, 

because of (jlj) and the forms of w and v. 

Using again the theorem of Moreau, it follows the conclusion of the theorem. 

□ 



5 Projection onto the monotone nonnegative cone 
via the projection onto the monotone cone 

There are various efficient methods for projecting onto the monotone cone emerging from 
the PAV algorithm (see e. g. PQ). These methods are intimately related to the special 
structure of the monotone cone, and their justification is more function theoretic than 
geometric. Due to their simplicity their usage is desirable for projection onto the monotone 
nonnegative cone too. 

As a consequence of Theorem [1] we have the following corollary; 

Corollary 1 Suppose that W is the monotone cone, that is, 

W = {x G R m : x 1 < x 2 < ... < x m }, 

and 

K = {x(ER m : < x 1 < x 2 < ... < x m } 
is the monotone nonnegative cone. Then for an arbitrary y G M. m it holds 

P K y = (P w y) + , 

where + stands for the lattice operation defined by the order induced by the nonnegative 
orthant in W m . 
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To exploit the efficiency of PAV-type algorithms in projecting onto the monotone 
nonnegative cone, we can proceed as follows: For an arbitrary y G M m we can determine 
PwV = (x 1 , ...,x m ), the projection of y on the monotone cone W, by using a PAV-type 
algorithm (e.g. the algorithm in [3]). Then, take the vector 

v = (x 1+ ,..,x m+ ), 

where x l+ denotes the "positive part" of the coordinate x l . Then, according the above 
corollary, the projection of y on the monotone nonnegative cone K C W engendering W 
is given by 

PkV = v. 
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